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Abstract 

For a general class of unitary quantum maps, whose underlying classical phase space is divided into 
several invariant domains of positive measure, we establish analogues of Weyl's law for the distribution 
of eigenphases. If the map has one ergodic component, and is periodic on the remaining domains, we 
prove the Schnirelman-Zclditch-Colin de Verdiere Theorem on the equidistribution of eigenfunctions 
with respect to the ergodic component of the classical map (quantum ergodicity). We apply our main 
theorems to quantised linked twist maps on the torus. In the Appendix, S. Zelditch connects these 
studies to some earlier results on 'pimpled spheres' in the setting of Riemannian manifolds. The 
common feature is a divided phase space with a periodic component. 

1. Introduction 

The quantization of an invertible map $ on a (i-dimensional compact manifold M. produces a unitary 
N x N matrix, the quantum map Un{^)- If there is a quantization recipe that works for an infinite 
sequence of integers N, one natural question is whether dynamical properties of the classical map $ 
are recovered in the semiclassical limit N — > oo. In this paper we prove that this is indeed possible for 
maps whose phase space M. is divided into several ^-invariant sets of positive measure, provided the 
quantization recipe satisfies a correspondence principle: quantum and classical evolution of observables 
must be equivalent in the semiclassical limit. The precise formulation of the necessary axioms is given 
in Section |21 These are kept fairly general and allow applications outside of quantum mechanics. In 
particular, A4 is not required to be symplectic. 

Let ifj S (j = 1, . . . , N) be an orthonormal basis of eigenvectors of Un(&), and 9j e M the 
corresponding eigenphases defined by the relation 

(1.1) U N ($) ifj = e 2 ^ tpj. 

In Section \5\ we prove an analog of Weyl's law |431 1441 [21 122] for the limiting distribution of the 
eigenphases as N — ► oo. Let us denote by /i an invariant probability measure of on M. which 
is absolutely continuous with respect to Lebesgue measure. We furthermore assume that the map is 
periodic on a collection of sets T>\, T>2 , . . . with positive measure [i{Pv) > and boundary whose fractal 
(i-dimensional Minkowski content is zero 1 , and that on the remaining set T>q the periodic orbits of <£> 
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Royal Society Research Grant 22355, and the EC Research Training Network (Mathematical Aspects of Quantum Chaos) 
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1 Fractal Minkowski contents and dimensions also play an important role in the Weyl-Berry conjecture for the distri- 
bution of eigenvalues of the Laplacian of domains in M. with fractal boundary, see 4 5 27 33 and references therein. 
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Figure 1. Two typical quantum eigenstates of a linked twist map on the torus for 
N = 201, one localised in the lower left quadrant Di (left image), the other in the com- 
plement D (right image). According to Theorem 17.11 approximately 25% resp. 75% 
of all eigenstates behave in this way — the latter states are in fact equidistributed on 
D , cf. Theorem IO 




FIGURE 2. Two untypical quantum eigenstates of the same linked twist map as in 
Figure ^ with N = 201. These eigenstates are localised near the boundary between 
the domains Do and Di; the eigenstate on the right has comparable mass in both 
domains. Eigenstates of this type form a sequence of density zero, cf. Theorem 17. II 
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Figure 3. The integrated density of quantum eigenphases In (9) of the linked twist 
map in Figured with N = 201. The jump at 9 = by \ reflects the fact that 25% of 
eigenphases are asymptotically zero. The remaining 75% of eigenphases are uniformly 
distributed (modulo one) in the interval [— i, i], in accordance with Weyl's law Q1.4|) . 
cf. Theorem 15 .21 

of any given period form a set with d-dimcnsional Minkowski content equal to zero. The number of 
eigenphases in some generic 2 interval [a, b] with b — a < 1 is then 

(1.2) lim = 1, . . . , JV : 0j G [a, b] mod 1} = V fi(V v ) / p v {6) dO 

N^ooN ^ J a 

where po(9) = 1 is the uniform probability density mod 1 and, for v > 1, p v (ff) is the uniform 
probability density supported on the points (k/n„) + ax> v ] here n v denotes the period of $ on 2?„, 
and the constant av v G K depends on the chosen quantization recipe. This means in particular that 
if A4 = T>o, the spectrum of the quantum map is uniformly distributed on the unit circle. A formula 
analogous to ljl.2|l has been obtained by Zelditch > Theorem 3.20, in the case of the wave group for 
a compact Riemannian manifold. A special case of formula i|1.2[) is proved in |49) for quantised contact 
transformations, and in |Sj for perturbed cat maps; in both cases the set of periodic orbits has measure 
zero and thus M. — T) . 

The main result of this paper describes the semiclassical distribution of the quantum map's eigen- 
states on the classical manifold. If the quantization constants ax> 1 ,Oix> 2 , ■ ■ ■ are linearly independent 
over the rationals, then, for N large, approximately N x p(D u ) of the N eigenstates localise on the 
set T) v (Section EJ. If in addition the classical map acts ergodically on T>q, then almost all of the 
eigenstates localised on T>o are in fact equidistributed on 2?o (Section [SJ. This latter result may be 
viewed as an extension of the Schnirelman- Zelditch- Colin de Verdiere Theorem, originally formulated 
for completely ergodic Hamiltonian flows (401 1451 1"HT1 1231 l2"Tl I50| and maps |%l 1481 14*51 115) . to maps with 
partially ergodic phase space. The possibility that eigenfunctions localise exclusively on the ergodic or 



Generic means here that a, b are not in the singular support of the limit density (which in the present case is a 
countable set). 
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on the integrable component had been conjectured by Percival in the 1970s [j""""j, and is known as the 
semiclassical wave function hypothesis. 

Two typical and two untypical examples of eigenstates of a quantised linked twist map on a two- 
dimensional torus (see Section """"2*}l . with N = 201, are displayed as a Husimi density plot 3 in Figures 
[""]and[""] The classical map acts as the identity on the lower left quadrant, and ergodically on the 
L-shaped complement fSection "3.3[) . Figure shows the integrated density 



Very similar quantum maps have been investigated numerically in |30| ; our theorems apply to those 
cases with sharply divided phase space and rational frequencies in all elliptic islands. A further inter- 
esting family of toral maps with mixed dynamics are the lazy baker maps |28j . Here the phase space 
is divided into countably many rational elliptic islands of total measure one; the hyperbolic dynamics 
takes place on a fractal set of Lebesgue measure zero. 

In the case of Hamiltonian flows with partially ergodic phase space 4 , Schubert [""""] has proved a result 
analogous to our quantum ergodicity theorem f Theorem 18. 1(1 . which however only holds for quasimode 
solutions of the Schrodinger equation and not necessarily for the eigenstates itself. The problem is 
that near-degeneracies in the spectrum of the quantum Hamiltonian (which in general cannot be ruled 
out) may create eigenstates that are extended across the entire phase space although the corresponding 
quasimodes are localised on the flow invariant components. For the quantum maps considered here, the 
spectrum can be controlled sufficiently well to rule out this possibility. A similar observation results 
from the analysis of the wave group on Ricmannian manifolds where geodesic flow is periodic on some 
open invariant component, see the Appendix by S. Zelditch. 

It should be emphasised that our results do not address the question of the possible localization of 
eigenstates on sets of measure zero (such as scarred eigenstates, bouncing ball modes or the recently 
discovered hierarchical states [21] )• Results in this direction have recently been obtained in the case of 
cat maps [""] Ell Ell and piecewise affine maps on the torus [T2j ; which neatly complement the proofs 
of quantum unique ergodicity for cat maps [El ESI EE] , parabolic maps """"Q and the modular surface 
gaEHl (see the survey [33]). 

Concrete examples of maps satisfying the axioms set out in SectionElare toral linked twist maps, see 
Sections 13. II and 1X2*1 The specific structure of these maps permits a more detailed asymptotic analysis 
which will be presented elsewhere j35j . 

The axioms in Section [51 are in fact sufficiently general to allow also applications to sequences of 
unitary matrices without quantum mechanical interpretation. In Scction "3.4l we discuss an application 
arising in the discretization of classical evolution operators. 



We thank A. Backer, M. van den Berg, E. Bogomolny, Y. Colin de Verdiere, S. De Bievre, F. 
Mezzadri, S. Nonnenmacher, J. Robbins, Z. Rudnick, R. Schubert, F. Vivaldi, A. Voros and S. Zelditch 
for stimulating discussions. 

2. Set-up 

Let M. be a <i-dimensional compact smooth manifold, and /i a probability measure on M which is 
absolutely continuous with respect to Lebesgue measure. 

3 See e.g. |34l IT) for a detailed discussion of Husimi functions. 

^Examples of such systems are the billiard flows in Bunimovich's mushrooms 1101 . 
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We fix an atlas of local charts <j>j : Vj — ► K d , where the open subsets Vj cover M.. In the following 
we thus identify subsets S of M. with subsets E of K d in the standard way. Let E be a subset of M. d , 
and 

(2.1) E(e) = UeM d :d(£,E)<e} 

its closed e-neighbourhood, where d( ■ , ■ ) is the euclidean metric on M. d . The s- dimensional upper 
Minkowski content of E is defined as 

(2.2) 9JT S (E) :=limsup(2e) s -'V(E(e)), 

where v is Lebesgue measure, see [22 Section 5.5]. We say E has Minkowski content zero if 9Jl* d (E) = 0. 
This is equivalent to saying that for every 8 > we can cover E with equi-radial euclidean balls of 
total measure less than 6. We say a subset S of M has Minkowski content zero if each of the sets 
Ej := 4>j(^\vj) C M rf has Minkowski content zero. 

We consider piecewise smooth invertible maps $ : M — > M. which preserve /i. By piecewise smooth 
we mean here and in the following that there is a partitioning of M. into countably many open sets Ui, 
i.e., M — \J i Ui and Ui HUj — 0, so that $|^. is smooth, 5 and the boundary set lj i dUi has Minkowski 
content zero. We will refer to this set as the domain of discontinuity of and call its complement 
M, — dUi the domain of continuity of $. 

Let Mjv(C) be the space of N x N matrices with complex coefficients. For a given infinite subset 
{index set) IcN, we say two sequences of matrices, 

(2-3) A := {yljvj-jveXi B := {-Bjv}jveii 

are semiclassically equivalent, if 

(2.4) \\A N -B N \\^0 

as N G X tends to infinity, where || • |j denotes the usual operator norm 

(2.5) \\A\\:= sup I'M. 

«/>ec«— {o} WWW 

We denote this equivalence relation by 

(2.6) A~B. 
Lemma 2.1. If A ~ B then Tr A N = TrB N + o(N). 
Proof. We have 

(2.7) ^\Tr A N -TrB N \ < \\A N - B N \\ -> 0. 

□ 

Let us define the product of two matrix sequences by AB — {Ajv-B/vj-jveZ; the inverse of A by 
A -1 = {Aj^j-jveii an< i its hermitian conjugate by A^ = {A' N }Nex- 

Axiom 2.1 (The correspondence principle for quantum observables). Fix a measure fi as above. For 
some index set IcN, there is a sequence Op := {Op n }nei of linear maps, 

Op N : C°°(X) -► Mjv(C), a >-> Pjv (a), 

so that 

(a) forallaeC°°(Al), 

Op(5)~Op( a )t; 

(b) for allai,a 2 € C°°(M), 

Op(ai)Op(a 2 ) ~ Op(aia 2 ); 



This means that and all its derivatives are bounded continuous functions Ui — * A4; we allow for the possibility 

that those bounds are not uniform in i. 
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(c) forallaeC°°(.M), 

lim — TrOp^fa) = / adii. 
N^ooN FNy J M 



N^oo N 

Examples of quantum observables satisfying these conditions are given in Section 13.11 In standard 
quantization recipes (such as the one discussed in Section [3.1 1) one in addition has the property that 

(2.8) Op( ai )Op(a 2 ) - Op(a 2 )Op( ai ) ~ Op({ai, a 2 }) 

where { , } is the Poisson bracket. This assumption is however not necessary for any of the results 
proved in this paper. The axioms (a)-(c) in fact apply to examples without quantum mechanical 
significance. One interesting case arises in the discretization of linked twist maps, where 

(2.9) Op(ai)Op(a 2 ) = Op(aia 2 ) = Op(a 2 )Op(ai), 
see Section l3~4l 

Axiom 2.2 (The correspondence principle for quantum maps). There is a sequence of unitary matrices 
[/($) := {Un($)}n£x such that for any a £ G°°{M) with compact support contained in the domain 
of continuity of 4>, we have 

J7($)- 1 Op(a)J7($) - Op(a o $). 

In the following we consider maps <& that may be periodic on a collection of disjoint sets T> v C A4 
[y = 1,2, . . .) of positive measure /z(2? M ) > 0, with periods n v so that $""1^ = id. In addition to 
Axioms [2.11 and we will here stipulate that there are constants ax>„ G K such that 

(2.10) [/($)""Op(a)~e(n„acjOp(a), e(z) := exp(27riz), 

for any a £ C°°(A4) with compact support contained in T> v and the domains of continuity of <&, $ 2 , . . . , 
Whereas Axioms l2~T1 and l2~2l are satisfied by all standard quantization schemes, condition (|2.10() is more 
restrictive: the constant old u could for instance be replaced by Op(/3„) where [3 V is a non-constant 
smooth function on T> v ] note the relation 

(2.11) [/($)"" Op(<z) ~ e(n v Op(/3 v )) Op(o) 

is still consistent with Axiom l2~2l Condition H2.1t)|l is however essential in the proofs of our main results 
since the spectrum of e(Op(/3„)) may, in general, be dense on the unit circle in C. 

In Sections 13.21 and 13.41 we discuss examples of semiclassical sequences of quantum maps satisfying 
the above Axioms 12.11 12.21 and condition (|2.1U|) . 

3. Example: Linked twist maps 

In this section we construct a well known example of quantum observables on the two-dimensional 
torus M. = T 2 := M 2 /Z 2 satisfying Axiom l2~Tl fcf. |16p. and corresponding examples of quantum linked 
twist maps satisfying Axiom l2~2l 

3.1. Quantum tori. It is convenient to represent a vector -0 £ C N as a function ip : Z/iVZ — > C. Let 
us define the translation operators 

(3.1) MKQH<MQ + i) 

and 

(3.2) MKQ) = e N (Q)ip(Q), 
where e]y(x) :— e(x/N) — exp(2irix / N) . One easily checks that 

(3.3) t^t™ 2 = t 2 ™ 2 ^ 1 ejv(TOim 2 ) Vmi,m 2 eZ. 

These relations are known as the Weyl-Heisenberg commutation relations. For m — (mi,m 2 ) £ Z 2 put 



(3.4) T N (m)=e N (^) 



,rn 2 j.mi 
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Then 

(3.5) T N {m)T N {n) = e N ( \ T N ( m + n ) 



with the symplectic form 

(3.6) u>(m, n) = m\n% — m%ri\. 
For any a £ C°°(T 2 ), we define the quantum observable 

(3.7) Op N (a) = a(™)T N {m) 

mez 2 

where 

(3.8) a(m)= [ a(£)e(-£ • m) d£ 

Jt 2 

are the Fourier coefficients of a. The observable Op N (a) is also called the Weyl quantization of a. 
Axiom im fa) is trivially satisfied. Axioms \'2. II (b) and (c) follow from the following lemmas. 

Lemma 3.1. For allai,a 2 £ C°°(T 2 ) 

(3.9) || Op A ,(ai)Op Ar (a 2 ) - Op N (a ia2 )\\ < —( ^ ||m|||a l( m)| ^ ll"ll|a 2 (" 

Proof. Using the commutation relations (|3.3|l we find 

(3.10) Pjv (ai)0 PAr (a 2 ) - ^ Si(m)a 2 (n)T N (m)T N (n) 



(3.11) = eN [ 2 ^(^^(m + ra) 

(3.12) = J2 e N (^^)d\( m )a 2 (k~ m )T N (k) 



m,k£Z, 2 

with fe = n + m, and hence 



(3.13) || Op JV (ai)Op A r(a 2 ) - Op Ar (aia 2 )|| < 
The lemma now follows from 

(3.14) \e(x) - 1| < |2ttx|, |w(m, n)| < ||m|| ||n 
Lemma 3.2. For any a € C°°(T 2 ) and R>\ 

(3.15) i-TrO Pjv (a)= / ad^ + O a . R {N- R ). 
Proof. Note that 



|al(m)l |(T 2 (n) 



(3.16) TrTjvM 



AT if m = 0mod A^Z 2 , 
otherwise. 



□ 



The lemma now follows from the rapid decay of the Fourier coefficients a(m) for ||m|| — > oo. □ 
Note that we have the alternative representation for Op N (a), 

(3.17) [0pNia ^ m = J2a(m,Q + ^ip(Q + m) 
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where 

(3.18) a(m,q)= / a(p 7 q) e(—pm) dp, 

JR/Z 

which is sometimes useful. In fact (|3.17|) permits to quantise observables a which are discontinuous in 
the g-variable. Note that if a is a smooth function of p and, for any v > 0, -j-^a(p, q) is a bounded 
function on T 2 , then, for any R > 1, there is a constant Cr such that 

(3.19) \a(m,q)\<C R (l + \m\)- R 

for all m, q. This fact is proved using integration by parts. Of course l|3.19|l holds in particular for 
smooth observables a € C°°(T 2 ). 

3.2. Quantum linked twist maps. A twist map ^> f is a map T 2 — » T 2 defined by 

(3.20) *' : @- > ( P+ / (9) ) m ° dZ2 

where / : R/Z — > R is piecewise smooth, i.e., the domain of discontinuity of / in R/Z has Minkowski 
content zero, cf. Section [21 Piecewise smooth functions of this type may be realised by taking a 
countable set of points = £o < £i < ■ ■ • < Coo = 1 with finitely many accumulation points in [0, 1], 
and assume that / G C°°([£j, for alH = 0, 1, . . . , oo. 

Obviously Lebesgue measure dfi = dp dq is invariant under f . A linked twist map $ is now obtained 
by combining two twist maps, ^ f t and *f?f 2 , by setting 

(3.21) $=Rof Jl oR- 1 o$ /j 
with the rotation 

(3.22) R:f P )^f q ) modZ 2 . 



Since M^, ^ f 2 and R preserve /z, so does $. More explicitly, we have 



(3.23) Ro^oR" 1 : J h P I mod Z 2 



and thus 

(3.24) <!- : ('') • ■ ( P A f l iq )r^) 



q) \q-h{p + f2(q)) / 
We define the quantization of the twist map 'J/ by the unitary operator 



(3.25) [U N (*fM(Q) = e 



-NV 



i>(Q) 



where V is an arbitrary choice of a piecewise smooth function R/Z — > R satisfying / = — V and 
V € C°° on the domain of continuity of /. 6 

Proposition 3.3. for any a € C°°(T 2 ) wit/i compact support contained in the domain of continuity 
of ^ f , we have 

(3.26) WUN^f)- 1 Op N (a)U N (* f ) -Op N (aoM, f )\\ = 0(N~ 2 ) 

where the implied constant depends on a. 



^In the case of smooth twist maps : T 2 — > T 2 it is more convenient to view / as a C°° function K/Z — > K/Z, 
thus avoiding the introduction of artificial discontinuities in /. The potential V may now be defined as a C°° function 
V : K — ► R with / = -V locally, and the condition that NV(Q + 1) = NV(^) mod 1 for every Q,N. Examples are 
f(q) = 2rq with t £ Z, for which V(q) = —rq 2 . In this case the correspondence principle stated in Proposition 13 , 3l is in 
fact exact, i.e., the right hand side of 13.261 is identically zero. 
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Proof. We have 

(3.27) [CM*/)" 1 Pjv (a)IM*/M(Q) 
and 

< — , / O m. \ 

mf 



N 



V 



Q + m 
N 



(3.28) 

since 

(3.29) 

Therefore 

(3.30) 

with 
(3.31) 



Op N (ao* f ) = 



Q m 
a ' m 'iV + 2iV 



(Q m 

I TV 27V 



1 N 



ip(Q + to), 



ip(Q + m), 



(a o */)(to, q) = e[mf(q)} a(m, q). 
HEM*/)" 1 Pjv (a)CM*/) - Op N (a o */)|| < max £ 



a to, g + 



2/V/ 



Cm (9, iV) = e {-N [V (q + ^) - V (q)} } - e [to/ (g + ^ 



Since |c m (q, AT)| < 2 and |a(m, g)| < (1 + |ra|) 5 , we have 



(3.32) 



max 



E 

|m|>JVV2 



" ( T^C^/V) 



< JV~ 



Let us denote by CS the projection of the compact support of a onto the q axis. CIS is a compact set 
which is in the domain of continuity of /. 

For |m| < TV 1 / 2 , Taylor expansion around x = q + yields (the second order terms cancel) 



(3.33) 
(3.34) 



V 



X+ 2n)- V ( X -2n)- V{x) N +0 {n^ 



( rr? 



N 



\N 3 



uniformly for all \m\ < N 1 / 2 and all q £ CS, provided TV is sufficiently large so that [q— N x l 2 , q+N x / 2 ] 
is contained in the domain of continuity. Hence in this case 



(3.35) 

and 

(3.36) 

(3.37) 



c m (q,N) = O 



( TO 3 

[n 2 



max 



E 



|m|<JVV2 



"'< m,q+ — )c m (q,N) 



< 0(Ar- 2 )max V 



771 &l m,q+ 2N 



= 0(N- 2 ). 



□ 



Proposition 3.4. Suppose V(q) = v = const for all q in some open interval I C R/Z. Then, for any 
T > 1 and any a G C°°(T 2 ) with compact support contained in R/Z x I C T 2 , we have 

(3.38) l|£M*/) Op N (a) - e(-Nv) Op N (a)\\ = 0(N- T ) 

where the implied constant depends on a and T. 
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Proof. We have 

(3.39) [{U N (* f ) - e(-Nv)) Op N (a)^](Q) 

'Q 



E~( Q . m 
a { m >N + 2N 



-NV 



N 



e(-Nv) \ip(Q + m). 



As before we split the sum into two terms corresponding to |m| < N 1 / 2 and \m\ > N 1 / 2 . For N large 
enough, the first term vanishes since a(m, q), as a function of q, is compactly supported inside the open 
interval /. The second term is bounded by 



(3.40) max V 



fry 

A m > 1 + ™ ) i e l- N V(Q)} ~ e(-Nv)} 



« a ,r N~ 



2N , 

for any T > 1 due to the rapid decay l|3.I9[l . □ 



The discrete Fourier transform Tn is a unitary operator defined by 

1 N ^ 

(3.41) [Fn4>]{P) = t=E ip(Q)e N (-QP). 



N ~=o 



Its inverse is given by the formula 

N-l 

(3.42) = "Ttt Y] 4>{P)e N {PQ)- 



N - ^ 

v JV p=o 

Proposition 3.5. For any a <E C°°(T 2 ) 

(3.43) Op N (a)F N = Op N (a o R) 
with the rotation R as in (I3.22|l . 

Proof. This follows from the identities F^tiJ-N = t% and T~^t 2 T^ — t\. □ 

The Fourier transform may therefore be viewed as a quantization of the rotation R which satisfies 
an exact correspondence principle, cf. Axiom |2~21 

The quantization of the linked twist map is now defined by 

(3.44) U N ($) =T N U N {^ h )^ 1 U N {^ h ). 

Proposition 3.6. For any a £ C°°(T 2 ) with compact support in the domain of continuity of <£>, we 
have 

(3.45) \\U N {<$>)- 1 0v N {a)U N {<S>)-0v N {ao<$>)\\=O{N- 2 ) 
where the implied constant depends on a. 

Proof. Apply Propositions 13.31 and 13.51 □ 

The quantum map £//v($) thus satisfies Axiom 12721 

Proposition 3.7. Suppose V\(p) = v\ = const for all p in some open interval I\ C R/Z, and V 2 (q) = 
i>2 = const for all q in some open interval I 2 C R/Z. Then, for any T > 1 and any a € C°°(T 2 ) with 
compact support contained in the rectangle I\ X I 2 C T 2 . we have 

(3.46) \\U N ($) Op N (a) - e[-N{ Vl + v 2 )\ Op N (a)\\ = 0(A^ T ) 
where the implied constant depends on a and T. 

Proof. Apply Propositions 13 . 41 and 13.51 □ 

In the examples considered in Section 13.31 we have v± — v 2 — 0, and furthermore, for any n =/= 0, 
acts as the identity precisely on the rectangle 1% x I 2 . Hence condition (|2.10() [cf. condition (c) of 
our central Theorem 15.21 below] is satisfied in this case. In cases where v\ + v 2 7^ mod 1 one may 
consider subsequences of N — > 00 for which {N(v\ + v 2 )} — > a, for any suitable fixed a € [0, 1]; here 
{ • } denotes the fractional part. 
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3.3. Ergodic properties of linked twist maps. The ergodic properties of linked twist maps are well 
understood |1 11 137j . Let [a»,6j] (i = 1,2) be subintervals of R/Z, and choose functions /, : R/Z — > R 
with 

(a) fi(q) = for q £ [a^bi], 

(b) fiifli) £ Z and fiibi) — fi(a>%) — fcj for some integer fc, 6 Z, 

(c) fi € C 2 ([a,i,bi\) with derivative /,-(<?) 7^ for all 9 G [aj,6j]- 

Let us define the constant 

(3.47) 7- i = sign(fc l ) max |/-(g)|. 

qe[ai,bi] 

Theorem 3.8. Suppose either of the following conditions is satisfied, 

(i) 7i72 < 0; 

(ii) Ifcil, \k 2 \ > 2 and 7l72 > C ~ 17.24445. 

Then the map l|3.21|l acts ergodically (with respect to Lebesgue measure fi) on the domain 

(3.48) V = {(p,q)eT 2 :pe [a u h]} U {(p, g) e T 2 : q e [a 2 ,b 2 ]}. 

The proofs of the two parts (i) and (ii) of this statement are due to Burton and Easton (Tl| , and 
Przytycki respectively. Both JI] and [SZj m fact establish the Bernoulli property for the action 
of $ on Vq under conditions (i), (ii). We expect that these properties hold under weaker conditions, 
e.g., for smaller values of Cq. The continuity of the map at the lines p = a\, b\ and q = a 2 , b 2 , assumed 
in condition (b), is probably also not necessary. 

For the linked twist map <& used in Figures we have chosen 

(3.49) [ 0l ,6i] = [02,62] = [il], 



(3.50) 

and 

(3.51) 



h(p) = 



-2p ifp€[|,l] 
if p ^ [| } 1], 



'2q ifge[|,i] 
if 



Vi(p) = 



V 2 (q) 



p 2 ifpe[i,l] 
ifp^[|,l], 



if ge [i 1] 
if^[|,l]- 



More explicitly, this particular map $ is obtained by first applying the twist 

if<ze [|,i] 



(3.52) 



followed by 



1 2" 
,0 1, 



p ' 



(3.53) 



'1 0' 

.0 1, 



'1 0^ 
.2 1, 



P 

v9y 



p ' 



if g ^ [|,1] 



ifpG [|,1] 



'1 (T 
.0 1, 



P 



ifp^Iil]- 

Clearly, this $ satisfies the above conditions (a)-(c) and (i). 
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3.4. Discretised linked twist maps. This section illustrates that the results of this paper may be 
applied to problems outside quantum mechanics, such as the discretization of maps on the torus, 
cf. [3] and references therein. A discretization of $ : T 2 — > T 2 is defined as the invertible 7 map 
$m : T 2 f — ► T 2 M with r f\ [ := (M _1 Z/Z) 2 , where $m is chosen in such a way that 

(3.54) lim sup d($ M (0,HO) =0; 

d( ■ , • ) denotes the Riemannian distance on T 2 . The discretised map induces a permutation matrix 
U N ($) :C N ->C N with N = M 2 defined by 

(3.55) [M*M(0 = ^ ° £ e T| f 

(we represent vectors ip £ C N as functions ip : T 2 f — » C). 

The "quantum" observables required in Axiom |2~TI are simply defined as multiplication operators, 

(3.56) [0 Pj¥ (o)V] (0=0(0^(0 
and therefore trivially satisfy Axioms [2. II fa) and (b). As to (c), 

(3-57) -TrO Pjv (a) = ^ £ a = / a d/i + O(M^). 

mi,m2- 



T 2 



It is easily checked that 

(3.58) Un^)- 1 Op N (a)U N ($) = Pw (ao$ M ) 
and thus 

(3.59) ||C/ A r($)- 1 Op Ar (a)C/ A r($)-Op Ar (ao$)|| = || Op N (a o $ M ) - PiV (a o $)|| 

(3.60) < sup |a($ M (0) - a(*(0)| 

(3.61) « sup d($M(0.*(0). 

and hence Axiom l2~2l is satisfied. Condition l|2.1U[) follows from a similar argument, with ap = 0. 

A concrete discretization of a linked twist map is obtained by replacing each twist map \& t by $ f M , 
with Jm{o) — M~ 1 \M f(q)], where [x] denotes the integer part of x. Note that R preserves T 2 M and 
therefore requires no further discretization. In this case 

(3.62) sup d($ M (0,$(0) =0(M~ 1 ). 

There is a simple geometric interpretation of the spectrum of Ujsr(<&). The map $m represents a 
permutation of N = M 2 elements, which can be written as a product of, say, v cycles C\ of length £j, 
i = 1, . . . , v. Each cycle corresponds to a periodic orbit of period li for the action of $a/ on T 2 f . Let 
£i be an arbitrary point on C%. An orthonormal basis of eigenvectors of Un{<&) is then given by the 
functions 

_ JiT^eUk/ii) if £ = $m (6) for some k = 0, . . . ,4 - 1, 
I otherwise, 

with eigenvalue Ay = e(J/£i), where i runs over the cycles and j over the integers 0, 1, . . . , l\ — 1. Hence 
in particular |</5y(£)| 2 = for every £ on the periodic orbit, and \<fij(0 \ 2 = otherwise. 

Due to the commutativity of the observables Op(a), the proofs of some of the statements in later 
sections may be simplified — especially those in Section 



^Invertability is not necessarily required in general discretization schemes. We assume it here to obtain a unitary 
representation. 
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4. Mollified characteristic functions 



Let us now return to the general framework of Section [21 Consider the characteristic function xv 
of a domain T> C M with boundary of Minkowski content zero. An e-mollified characteristic function 
Xv € C°°{M.) has values in [0,1] and xt>{x) = Xv(%) on a set of x of measure 1 — e. Since V 
has boundary of Minkowski content zero, we can construct such a smoothed function for any e > 0. 
Furthermore we are able to construct e-mollified xv whose support is either contained in 2?, or whose 
support contains T>, again for any e > 0. Note that if xv is e-mollified, so is Xv f° r any n £ N with 
the same e. 

After mollification, we may associate with a characteristic function xt> a quantum observable 
Op N (xv)- Since Op N (xv) is in general not hermitian, it is sometimes more convenient to consider the 
symmetrised version, the positive definite hermitian matrix 



The following proposition describes the distribution of eigenvalues of Op^ m (xx>), and suggests that 
the operator may be viewed as an approximate projection operator onto a subspace of dimension 
~ N x fi(V). 

Consider a sequence J := { Jn}ngi of sets Jn C {1, . . . , N}. The quantity 



provided the limit exists, is called the density of J. 

Proposition 4.1. Suppose xv is an e-mollified characteristic function, and suppose fij > (J = 
l,...,N) are the eigenvalues of Op^ m (xx>)- Then there are set sequences J := {Jn}nex and J' := 
{J' N }Nei with densities 




(4.3) 



A (J) := lim 




(4.4) A(J)= M p)+0(e 1 / 3 ), A(j') = l-M(^) + 0(e 1 / 3 ), 



such that 



(i) ft = 1 + Ofc 1 ^) for all j e J N ; 




(4.8) 




3=1 



Therefore 



(4.9) 




and thus 



(4.10) 




14 



JENS MARKLOF AND STEPHEN O'KEEFE 



( 4 - n ) *^E^ = ^ 

where Hm = {j : > 1/2}. By Chebyshev's inequality, i|4.11[) implies that 

(4.12) ^ £ i?Af : > 7} - 0(e/ 7 ), 
for any 7 > 0. This yields the bound 

(4.13) ton 1 £ W =0(7 1/2 + £ /7), 
since < /z/ < 1/2. So 

(4.14) lim 1 y 1= hm 1 V [(^ - l) 2 - + 2/x 



(4.15) =^jl E M + 2^]+0( e ) 

1 w 

(4.16) = lim - £ M + + 0(e) + 0( 7 1/2 + e/ 7 ) 



(4.17) =M^) + 0(e 1 / 3 ), 

if we choose 7 = e 2 / 3 , and hence the corresponding set sequence H :— {Hn}n<zj has density A(ff) = 
fi{V) +0(e 1 / 3 ). Once more in view of Chebyshev's inequality, (|4.1U|) implies 

(4.18) lim 6 Hat : - I) 2 > 5} = 0(e/«5). 

Choosing <5 = e 2 / 3 , this means that for a subsequence of j £ Hn of density /^(£>) + 0(e 1 ^ 3 ) we have 
/ij = 1 + 0(e 1 ^ 3 ). The corresponding result for j follows by the same argument from (|4.12() . □ 

5. Trace asymptotics and Weyl's law 
The following proposition is the key tool to understand the distribution of eigenvalues of [/at($). 

Proposition 5.1 (Trace asymptotics). Suppose $ : M. — > M. is piecewise smooth, and 

(a) 3> n |.p, = id on some set D' cM with boundary of Minkowski content zero; 

(b) the fixed points of $ n on T> := A4 — T>' form a set of Minkowski content zero; 

(c) there is a constant a-p' G R such that for any a e C oa (A4) with compact support contained in 
D' and the domains of continuity 0/$™, we have 

[/($)™Op(a) - e(nax>0 Op(a). 

Then 

(5.1) lim — TrU N ($) n = e(na v ,) n(D'). 

N^oc N 

Proof. Given any e > 0, we can find an integer R and a partition of unity on M by e- mollified 
characteristic functions, 

R 

(5.2) l = Xbad(£)+Xi>'(£)+E^) y Z £M 

r=l 

with the properties 

(i) the interior of the support of Xbad contains the domains of discontinuity of in M. for all 
v = 1, . . . , n, and all fixed points of in T>, and is chosen small enough so that J XbuddfJ. < e; 
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(ii) the support of \V' is contained in V and the domains of continuity of $ w for all v = 1, ... , n, 
so that fi(T>') ~ J xvdfi < e and all fixed points of <I> n are contained in the interior of the set 

SUpp Xbad U SUpp XV ; 

(iii) the support of Xr, with r = 1, . . . ,R, is chosen small enough, so that suppxrn3>™(suppxr) = 
for aH£eM. 

Properties (i) and (ii) are possible since the fixed points in T> and the discontinuities form sets of 
Minkowski content zero. To achieve (iii) note that the closure of K, = M. — (suppxbad U supp%i>/) 
does not contain any fixed points, and $ is continuous on /C. Hence there is a sufficiently small radius 
T] = 77(e) such that for all balls B v C K. we have B v H $ n (B v ) = 0. 
By the linearity of Op, we have 

R 

(5.3) TrU N (*) n = Tr[UN@) n Op N (xtod)] +Tr[U N ($) n Op N (xv>)] + £>[M$) n Pjv (Xr)]. 

r=l 

We begin with the first term on the right hand side: 

(5.4) Tr[U N ($) n CWxbad)] = Tr[U N ($) n Op^tXbad)] + o £ (N), 

with the symmetrised Op^ m (xbad) as defined in 14.1(1 . Suppose tpj and \ij > are the (normalised) 
eigenstates and eigenvalues of Op^ m (xbad)- Then 



N 



N N 

(5.5) |Tr[^($)"Op^ m (xbad)]| = ' 



3=1 3=1 

= TrOp s ^ m (xbad) = TrO Pw (xbad) + o 6 (iV) = 7VO(e) + o 6 (JV). 

Since by condition (c), U(^) n Op(xv) ~ e(naii')Op(xi)')i we find for the second term on the right 
hand side of (15. ! 



(5.6) Op^xuO] =e(na^)TrOp JV (x 2 ,0 + o,(iV) 

(5.7) =jVe(no!T.') / Xvdn + o e (N) 

(5.8) = 7V e ( nap ,) {/i(D') + 0(e)} + o e (JV), 

where we have used (ii) in the last step. 

For the last term in the sum (15.311 we have 



(5.9) E7(*) w Op(xr) ~ J7(*)"Op(xy 2 )Op(xy 2 ) ~ Op^ 2 o $-")[/(<D)"Op(xy 2 ) 
so 

(5.10) Tr^jv^)" PAr (xr)] - Tr[0 Pjv (xy 2 o $-")[/ JV ($)" Op^ 2 )] + o e (N) 
(5.H) = Tr[0 Pjv (xy 2 ) Op N (xl /2 o $-™)C7^($)«] + o e (7V) 

(5.12) = Tr[0 Pjv (xy 2 • x^ /2 ° $"")f/jv($)"] + o e (iV) 

(5.13) =o £ (iV) 

since Xr^ 2 ■ Xr^ 2 ° = in view of (iii). Therefore 

(5.14) TrU N ($) n = Ne(nav>){KV')+0(e)} + o e (N), 
i.e., 

(5.15) lim ±-TrU N (<S>) n = e{na v >) fi(V) + 0(e), 

N^oo iv 

which holds for every arbitrarily small e > 0. This concludes the proof. □ 
Theorem 5.2 (Weyl's law). Suppose $ : M. — > M. is piecewise smooth, and 
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(a) there are disjoint sets T> v (u = 1, 2, 3, . . .) with boundary of Minkowski content zero, on which 
$ is periodic with period n v > 1, i.e., 

and, for every < \n\ < n v , the fixed points of $" on D v form a set of Minkowski content 
zero; 

(b) for every n ^ 0, the fixed points of $™ on Dq = M. — U^Li T^v form a set of Minkowski content 
zero; 

(c) there is a constant ax> v £ ffi such that for any a e G°° {M.) with compact support contained in 
T> v and the domains of continuity of $, <& 2 , . . . , , we have 

E/(*) B "Op(a) ~ e{n v a-D v ) Op(o). 
Then, for every continuous function h : S 1 = M/Z — > C, 

-. N oo . 

(5.16) lim - E K0j) = E / d0 
where 

(5.17) p o (0) = 1, 



and, for v > 1, 
(5.18) 



1 { h \ 

■M = — E M e «^ • 



Here 

(5.19) «5 s i(0) = E <5(0 + m) 

mGZ 

denotes the periodicised Dirac distribution. 

Proof. For every v such that n„ divides n we have $ ra |p — id. A simple modification of the proof of 
Proposition 15 . II yields therefore, for ra^O, 

(5.20) lim — TiU N ($) n = V e(natD u ) Hp"). 

iV— too rv — ' 

n v \ n 

The only difference in the proof is that T>' is divided into the domains X?i,X?2, ■ ■ • with different inte- 
gration constants ax>„ ■ For every e > there is a K = K c such that 



(5.21) J (J Z>„ 



i/=K+l 



< e. 



Hence one effectively deals with only finitely many domains T>\,... ,T>k and shows, following the steps 
in the proof of Proposition l5.ll that 

1 K 

(5.22) lim — Tr U N (Q) n = V e(nav v ) /i(2? v ) + 0(e), 

AT— too iV 

!/=l 

n„ |n 

which in turn yields (|5.2U|) . 

Let us first assume that the test function h has only finitely many non-zero Fourier coefficients, i.e., 



(5.23) 



h{6) = J2Hn)e(n9) 
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is a finite sum. We then have 

N 



(5.24) Jim j7^2h(9j)= lim i ^ fc(n) Tr 

N^oo iv ^ — ' iv — »oo IV ' — ' 

oo 

(5.25) = £(0) + J] fe(n) ^ e(na©J p(V v ) 

n€Z-{Q} v=l 
n v \n 

oo 

(5.26) =h(0) + ^2fi{V u ) 53 Kn)e{na Vu ). 

v=l neZ-{0} 

Since 

(5.27) ft(n)e(napj = — ^ ^ W e KJ e ( — ) = — 53 h \ 



n u — 1 / ; \ i n ^ ~~ 1 

, s , / AC 

n^Z A:— n£Z x ' fc— 



we obtain 

( 5 - 28 ) J™ 4 E w = \ mPo)^(o) + 53 ^ 53 /> 



■-- ... . , ■ , . - , > h np 



which proves the theorem for /i with finite Fourier series. We now extend this result to test functions 
he C 1 ^ 1 ). Let 

(5.29) h K {9) = 53 h{n)e{n6) 

riGZ 
\n\<K 

be the truncated Fourier series. Since h e C 1 ^ 1 ), its Fourier series converges absolutely and uniformly 
and hence, for any e > 0, there is a K such that hx{9) — e < h{6) < hx{9) + e for all 9 e S 1 . By l|5.28l) . 
the limits of the left and right hand side of 

N N N 

(5-30) - 5] h K (6j) -e<^53 < ^ 53 M*i) + e 

exist and differ by less than 2e, hence (|5.28|l holds also for the current h. The extension of l|5.28|l to h in 
C(§ 1 ) is achieved by the same argument, i.e., by approximating h pointwise by functions h t e C 1 (S 1 ) 
sothatfc e (0)-e</i(0)<ft e (0)+e. □ 

6. Generalised Weyl's law 

Proposition 6.1 (Generalised trace asymptotics). Choose 3> and Un{<&) as in Provosition \5.1\ Then 
for every a £ C oa (M) and n 0. 

(6.1) lim ±- T Tr[Op N (a)U N ($) n ] = e(nav>) ( a dp. 

N^co iv Jx>i 

Proof. By linearity of the relation Q6.1[l we may assume without loss of generality that a is real and 
min^ a(£) > 0. This implies that a 1 / 2 € C°°(.M). Analogously to the proof of Proposition l5.il we have 

(6.2) Tr[0 Pjv (a)Z7jv($r] = Tr\U N ($) n Op w (x bad • a)] 

R 

+ Tr[U N ($) n Op N (xv • a)] + 53 Tr\U N ($) n Op N (xr ■ a)} + o t (N). 

r=l 

The proof is concluded in the same way as the proof of Proposition l5.il with all mollified characteristic 
functions \ replaced by x • a. □ 
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Theorem 6.2 (Generalised Weyl's law). Choose $ and Un{$) as in Theorem 15. gl Let tpj £ C N 

(j = 1, . . . , N) be an orthonormal basis of eigenstates o/t//v (<!>), with corresponding eigenphases 9j G § . 
Then, for every a £ C°°(A4) and every continuous function h : S 1 —* C, 



(6.3) 



J im vEM^)<Op JV (a)Vi > Vi>=X)/ «^ / h{p)p v {O)d0. 



Proof. We may assume again without loss of generality that a is real and min^ a(£) > 0. In view of 
Proposition 16. II and the proof of Theorem 15. 21 we have for every hpc with finite Fourier expansion (as 
in (E23) 



(6.4) 



^ N oo . 

J im My2 h K(&j)(Op N (a)<fj,ipj)=y] / 

7 = 1 l/=0 JI3 - 



ad^ / h K (0)p v (6)d6. 



For any ft, > we have 
(6.5) 



N N 

J2 KWOpjfWw, Vj ) - J2 h (»i)\\ Op A r(a 1 / 2 )^| 
i=i i=i 



< sup h 



Tr[0 Pjv (a) - Op JV (a 1 / 2 ) Op^a 1 ' 2 )*] = o(N) sup ft. 



Hence (16 .4|) is equivalent to 

-i N oo „ -i 

(6.6) w lim o -^ftx(^)||Op JV (a 1 / 2 V J || 2 =^ / adp h K {6) p v {9) d6. 



We now use the same approximation argument as in the proof of Theorem^ for ft e C 1 ^ 1 ). Given 
any e, there is a if such that hx(0) — e < ft(0) < hxiO) + e for all € S 1 . The limits of the left and 
right hand side of 



(6.7) 4 - £ ]H Op^Ca 1 / 2 )^!! 2 < l^ft(^.)|| Op JV (a 1 / 

3=1 J'=l 

(6.8) < -^[ft^(^) + e ]||0 Pw (a 1 / 2 )^|| 2 

i=i 

differ by less than 

1 W 1 

(6.9) 2esupft if lim — V || Op^a 1 / 2 )^ 2 < 2esupft,f lim — Tr[Op N (a^ 2 ) Op w (a 1 / 2 ) t ] 



i=i 



(6.10) 



2e sup ft#- / a dp 
Jm 



which can be arbitrarily small for e — ► 0. Thus 



(6.11) 



N oo - -1 

^EM^)iiop^(o i/2 )^ir=x; / ^ / mp»{Q)do. 



A similar approximation argument shows that (|6.11|l holds also for all continuous ft. In view of l|6.5f) . 
the relation (|6.11|) is equivalent to 1|6.3[) . The assumption ft > can be removed by using the linearity 
of (EH in ft. □ 
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7. Localization 



The sequence := {iPn}n£I of vectors ipN € C — {0} is said to be semiclassically localised in the 



domain T> if for every a £ C°°(A4) with a\x> = 0, we have 
(7.1) 



lim (gPjvWl^ = o. 

AT-i-oo W'J'W 



Theorem 7.1. Choose $ and Un(<&) as in Theorem \5.<l and let (fx, ... , <pn £ C N be an orthonormal 
basis of eigenstates ofUjf(<p)- Then there are set sequences J := {Jn}n&x and J' '■= {J'n}n£I with 
densities 

(7.2) A(J) = M (X>„), A(J') = l-M(^o), 

such that 

(i) for any a £ C°°(A4) with a\-p Q = 0, 

J im 4 E K°pjv( a ) ( ft> < ft)l 2 = °; 



(ii) /or any a S C°°(.M) a|jvi-x> = 0; 

J™ at 51 K°Pjv(a)<y3j,^)| 2 = 0. 



Proof. For any given e > 0, there is a constant A' = K e such that 



(7.3) 

Consider the following subset of S 1 = 

K n„-l 

(7.4) 6 £ := |J |J 



u ^ 



< e. 



/c e k 
h «-d„ - — — : 1- aD„ 



A"^ 



Z. 



i/=l fe=0 

We may construct a continuous function h — h e with values in [0, 1] such that 



(7.5) 



fc(0) = 



if 6» e e e 



Consider the set Ja^ £ of j, for which 0j ^ 62c- The corresponding set sequence J c = {</Ar, £ }Arez has, 
in view of Weyl's law (Theorem I5.2fl . density A(J £ ) = /i(£>o) + 0(e), where the implied constant is 
independent of A~, since the measure of G>2 £ is at most 

K n v -l 



(7.6) 
Now, 
(7.7) 

(7.8) 



EE 

v=l k=0 



4e 
Kn h 



4e. 



E l(°Pw( a )^"'^)| 2 ^ J2 h (0j) \{Op N (a)<pj,<pj}\' 



3=1 
iV 



<^M^)l|Op Ar (aV J || 2 

iV 

= X>(0i)(Op„(M> i ,^) + Oe (iV). 



(7.9) 

Hence, by Theorem 16. 21 

(7.10) limsup- Y, \(Op N (a)^,^)\ 2 <Y \a\ 2 dft / h(0) Pv (0)M. 



2(1 
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Now, under assumption (i) of the theorem, 



/ \a\ 2 dp = 0, 



(7.11) 

and furthermore 
(7.12) 

since h is supported outside the support of p%, . . . , pk- For the remaining sum, 



V f \a\ 2 dp [ h(9) p v (9) d9 = 



(7.13) 

by (|7~3) . Hence 
(7.14) 



^ / |a| 2 d/^ / h{6)p v {9)d9 < V / |a| 2 d^ 
limsup-^ ^ KOp^a)^-,^-)! 2 = 0(e) 



< emax lal 2 



with e > arbitrarily small. Therefore there is a sequence of values e = e_/v such that e^v — > as 
N — > oo, and 



(7.15) 



limsupi V |(0p A ,(a)^,^-)| 2 = 



where Jjv := Jn^ n . The proof for case (i) is complete. 

As to case (ii), define J' N e as the set of j, for which 9j G 9 e . The corresponding set sequence 
= {J' Ne } NeI has density A(J^) = 1 - p(V a ) + 0(e), recall Weyl's law ( Theorem ICT. Then, 
analogous to l|7.10|) . 



(7.16) 



limsup- V |<0 Pjv ( a )^,^)| 2 <^ / |a| 2 d M / [1-M0)]P,W^, 

AT— >-oo JV .~T7 .._ n JT>„ Jo 



where all but the v — term vanish, since a|»--D = 0. So 



(7.17) limsup — V |(0 PAr (a)^,^-)| 2 < / |a| 2 d^ / 



\a\ 2 dp / [l-h(9)]d9 = 0(e). 



Corollary 7.2. There are set sequences I := {In}n£X and I' := {I' n }n£X with densities 

(7.18) A(I) = p J (V ), A(l') = l-p(V Q ), 
such that 

(i) the eigenstates <pj , j G In, are semiclassically localised in T>q ; 

(ii) the eigenstates (fj , j G I' N , are semiclassically localised in M. — T>o . 

Proof. This is a straightforward consequence of Theorem 17. II and Chebyshev's inequality. 

Theorem 7.3. If in addition to the assumptions of Theorem 1 7. 1\ the phases a v (v — 1, 2, 3, . . .) 

linearly independent over Q, then there are set sequences J v := {J n }ngi with densities 

(7.19) A(J 1 ') = pb(T> v ) 
such that for any a G C°°(M) with a\x> v = 0, 

(7.20) lim 1 Y, \(Op N (a)^,^)\ 2 = 0. 



□ 



□ 
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Proof. This follows from a slight modification of the proof of Thcorcm l7.ll where the set Jfj is approx- 
imated by the set Jj^ e comprising those j for which 



(7.21) 9jG |J 



' r k eS k eS 



h o.j> , r o/.v v H 



mod 1. 



fc=0 

The crucial observation is that, by the linear independence of the a v over Q, there exists a S = 5 €> k > 
small enough, such that sets corresponding to different v = 1, 2, 3, . . . , K are disjoint. □ 

Thus, by Chebyshcv's inequality, for a subsequence of j G Jj^ of density [i(V u ) the eigenstates (fij 
are semiclassically localised in T) v . 

8. Quantum ergodicity 

Let us now turn to the question of quantum ergodicity for maps which have one ergodic component 
2?o and are periodic on a remaining countable collection 2?i,2? 2 , ... of domains. Examples of linked 
twist maps with this property are discussed in Section 13.31 

Theorem 8.1. Choose $ and £/jv($) as in Theorem \5.''A and suppose $ acts ergodically on T>q. Let 
cpx, . . . , ipN S be an orthonormal basis of eigenstates ofUnf(<p). Then, for any a g C°°(Ai), 



(8.1) lim 4 Y 



N~>co N 



{Op N (a)tpj,(pj) ad/j, 



2 



0. 



Proof. We may assume without loss of generality that J v ad/i = and \a\ < 1. It is then sufficient to 
show 8 

(8.2) S 2 (a,N) := 1 ]T \(Op N (a)^, ^)| 2 - 

as TV — > oo. 

For any given T > 1, we may write 

(8.3) a = a T + a' T + a^ ad 
where 

(i) ar € C°° has compact support contained in 2?o and the domain of continuity of $, $ 2 , . . . , $ T , 
and furthermore J aT d\i = 0, |ar| < 1; 

(ii) a' T € C°° is supported inside .M — Vq, and |a^| < 1; 

(iii) J M \a^\ 2 d^<T-\ 
By the triangle inequality, 

(8.4) S 2 {a 1 N) 1 ' 2 < S 2 (a T ,N) 1 / 2 + S 2 (a' T ,N) 1 / 2 + S 2 (a^ d ,N) 1 / 2 
By Theorem Ol 

(8.5) lim S 2 (a' T ,N) = 0. 
Furthermore, by the Cauchy-Schwartz inequality, 

1 N 

(8.6) S 2 (4 ad , AO < ^EH°Piv(4 ad )^|| 

i=i 

(8.7) = l T rO PAr (|4 ad | 2 )+ OT (l) 
and hence 

(8.8) limsu P S' 2 (a^ ad ,iV) < T~ l . 



8 The argument presented here is inspired by the proof of quantum ergodicity for cat maps 1141 1381 . cf. also 
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As to the remaining term, 

(8.9) S 2 (a T ,N) = 1 £ |(0 PAr (a T )^,^)| 2 , 

it remains to be proved that the limsup of l|8.9|) can be made arbitrarily small for sufficiently large T. 
To this end define the ergodic average of ar by 

T 

(8.10) 



Since (pj are the eigenfunctions of [//v($) we have 

T 

i = — , i , 

■11) 



(8.12) 
(8.13) 



(8.14) 
by Axiom Now 
(8.15) 

(8.16) 

(8.17) 
We have 
(8.18) 
where 
(8.19) 



S 2 (a T , N) = 1 £ ± £ <C^($)- W Op JV (a T )C/ JV ($)" W , ^> 

n=l 

j=l n=l 
JV T 

- a?E M*) _B Op w (a T )f/iv($)"^ 

3=1 "=1 
f ^ 2 

= wE II°Pjv (a?)^' || +ot(1), 



3=1 



f N 

JjYl ll°Pw( a T)^| 



2 1 



TV 



3=1 



3=1 
1 ^ 

3=1 

/ \a^\ 2 d^ + o T {l). 



I \al\ 2 dv) <([ \a T \ 2 d^) +([ \a h T MT \ 2 dv 

JV J \JVa J \Jv Q 

T T 

T 1 \ " „ x,n badT 1 \ " bad „ x.n 



1/2 



We have by Cauchy-Schwartz 



(8.20) 

(8.21) 
(8.22) 



badT 1 2 



T 

d /i= J_V I (4 ad o$^)(a^o$*)^ 

T 1 /2 

5 (X l4ad ° ** |a (X l4ad ° $fe|2 dM 

- { |4 ad | 2 d M <i 

■/Co J 



1/2 



using the $-invariance of /j, and assumption (iii). Therefore 

(8.23) limsup 5a (a T , JV) < / |a T | 2 G?/i + 0(T -1 ^ 
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Since $ acts ergodically on Vq, we have a mean ergodic theorem for test functions a € L 2 (.M), i.e., 





and hence lim sup jy^oo Siip-Ti N) becomes arbitrarily small for T sufficiently large. 
Corollary 8.2. There is a set sequence I := {In}nei with density A(J) = /i(£>o) such that 
(8.25) (Op N (a)ipj,ipj) -> / adfi 



□ 



□ 



Appendix A. Converse quantum ergodicity. By Steve Zelditch 



The purpose of this appendix is to briefly make some connections between the foregoing article of 
Marklof and O'Keefe and some results in our articles |4fil 147) concerning the distribution of eigenfunc- 
tions and eigenvalues on certain Riemannian manifolds with an open invariant component where the 
geodesic flow is periodic. 

At the same time, we wish to emphasize the relevance of all of the results to the converse quantum 
ergodicity problem: are quantum ergodic systems necessarily classically ergodic? As discussed in the 
preceding article, it is well-known that classical ergodicity implies quantum ergodicity. However, there 
are few results in the converse direction. There could exist non-ergodic classical systems with quantum 
ergodic quantizations, because the invariant sets might not be quantizable in any suitable sense. 

When the system is periodic in an open invariant set, one expects to be able to quantize this feature 
and prove that the system cannot be quantum ergodic. Examples where this has been carried out 
are given in the preceding article and in j32] (see also jlj)- Aside from fully integrable systems, these 
appear to be the only examples where quantizations of classically non-ergodic systems have been proved 
to be non-quantum ergodic. It is very plausible that many other (if not all) classically non-ergodic 
systems are quantum non-ergodic, e.g. quantizations of KAM systems, but no rigorous proofs of this 
exist at this time except in the partially periodic case. 

The examples studied in |47| were special Riemannian manifolds which we might informally call 
'pimpled spheres': We take the standard sphere (S n ,go) and deform the metric, and possibly the 
topology, in a polar cap B r (xo) of some small radius r to obtain a new Riemannian manifold (M,g). 
We assume that M r := M\B r (xo) = S r l := S n \B r (xo) as Riemannian manifolds. We then denote 
by A4 r C S*M the smooth manifold with boundary formed by the closed geodesies of (S n ,go) which 
lie in M r . The boundary consists of closed geodesies which intersect dB r (xo) tangentially. Thus, the 
boundary between the invariant periodic component and the remaining component is smooth and in 
fact the geodesic flow is 'almost-clean' in the sense of jJSj. This condition resembles (but is more 
restrictive than) the Minkowski content zero condition on the boundaries of the mixed phase space 
components in the preceding paper. Also, we did not fix the metric or the type of dynamics of the 
geodesic flow in the non-periodic component. In 0Bj Theorem 3.20, we used the method of moments 
to determine the Szego limit measure of the wave group in these examples, i.e. the limit measure given 
in Marklof-O'Keefe's Theorem 15. 21 As in the their paper, we used these results to obtain eigenfunction 
distribution results. The following is a corollary of Theorem A (b) of proved in [47) : 

Proposition A.l. The Laplacian A of a pimpled sphere (M,g) is never quantum ergodic. 

Let us sketch the proof combining the argument of [47) and that of the foregoing paper. We assume 
the reader is familiar with pseudodifferential and Fourier integral operators, which are the observables 
and quantum maps in the Riemannian setting. 

Proof. Let A be a zeroth order pseudodifferential operator on M with essential support in M r , i.e. 
assume the complete symbol of A is supported in M. r . Since the geodesic flow of (M,g) is periodic of 
period 2tt in A4 r , Ae 2m ^ is a zeroth order pseudodifferential operator on M. Using the calculation 
of the symbol of the wave group in 18 , it is simple to see that the principal symbol of Ae 2my ^ equals 
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e la UA for a constant e IQ (a Maslov phase). Here, a a is the principal symbol of A. Thus, e 27rl ^ satisfies 
the hypothesis 12.10fl of the preceding article. 

Now argue by contradiction. If A were quantum ergodic, we would have 



(A.l) (Bfy^j) / <r B d£ 

J S'M 

along a subsequence of density one of any orthonormal basis {<fij} of A-eigenfunctions. Here, C is 
normalized Liouville measure and B is a zeroth order pseudodifferential operator. Letting B = Ae 27T1 ^~^ , 
we would obtain 

(A.2) (Ae 27riVS 0j,^-) ^e la / a A dC. 

However, we also have 



S'M 



(A.3) (Ae 2 ^^ , = e 2 ^ (Afo , ~ e 27riA ^ / a A dC. 

J S'M 

Here, {Aj} are the eigenvalues of y/~K associated to (f>j. It follows that e 27rlAj — * e la along a subsequence 
of eigenvalues of density one. 

But this implies that all of the eigenvalues of e 27nv/ ^ cluster around e 1Q , i.e. that 

(A.4) dfi X := jyTvT E ^ 

Here, N(X) — #{j : Xj < A}. But the weak limit of d(i\ was calculated in Theorem 3.20 (cf. also 
Theorem 15 .21 in the preceding article), and shown to be 

CiMrl ( C(Mr) \ ^ 

(A - 5) £^* ete H £(^WJ 

This contradiction proves that (M, g) is not quantum ergodic. 

□ 

In comparison, the argument in |47) used a trace formula for TrAH\ where is the orthogonal 
projection onto the span of (f>j with Xj < A. The above argument based on individual elements seems 
more vivid. 
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